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(57) Abstract: A method, system, and compater program product related to the diagnosis of diabetes, and is directed to predicting 
the long-tenn risk of hyperglycemia, and the long-term and short-term risks of severe hypoglycemia in diabetics, based on blood 
glucose readings collected by a self-monitoring blood glucose device. The method, system, and computer program product pertain 
directly to the enhancement of existing home blood glucose monitoring devices, by introducing an intelligent data interpretation 
component capable of predicting both HbAic and periods of increased risk of hypoglycemia, and to the enhancement of emerging 
continuous monitoring devices by the same features. With diese predictions the diabetic can take steps to prevent the adverse conse- 
quences associated with hyperglycemia and hypoglycemia. 
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METHOD, SYSTEM, AND COMPUTER PROGRAM PRODUCT 
FOR THE EVALUATION OF GLYCEMIC CONTROL IN 
DIABETES FROM SELF-MONITORING DATA 

5 

CROSS-REFERENCES TO RELATED APPLICATIONS 

The present invention claims priority from U.S. Provisional Patent Application 
Serial No. 60/193,037 filed March 29, 2000, entitled "Algorithm for the Evaluation of 
10 Glycenuc Control in Diabetes From Self-Monitoring Data" the entire disclosure of 
which is hereby incorporated by reference herein. 

US GOVERNMENT RIGHTS 

This invention was made with United States Government support under Grant 
15 Nos. NIH / NIDDK: ROl DK 28288 and NIH / NIDDK: ROl DK 5 1 562, both 

awarded by National Institutes of Health. The United States Go vemment has certain 
rights in the invention. 

FIELD OF THE INVENTION 
20 The present system relates generally to Glycemic Control of individuals with 

diabetes, and more particularly to a computer-based system and method for evaluation 
of predicting glycosylated hemoglobin (HbAu and HbAi) and risk of incurring 
hypoglycemia. 

25 BACKGROUND Qf THE INVENTION 

Extensive studies, including the Diabetes Control and Complications Trial 
(DCCT) (See DCCT Research Gio\q>:*The Effect Of Intensive Treatment Of Diabetes 
On The Development And Progression Of Long-Term Complications Of Insulin- 
Depcndent Diabetes Mellitus. New England Journal ofMediciney 329: 978-986, 

30 1 993), the Stockholm Diabetes Intervention Study {See Reichard P, Phil M: Mortality 
and Treatment Side Effects During Long-term Intensified Conventional Insulin 
Treatment in the Stockhohn Diabetes Intervention Study. Diabetes, 43 : 3 1 3-3 1 7, 
1994), and the United Kingdom Prospective Diabetes Study {See UK Prospective 
Diabetes Study Group: Effect of Intensive Blood Glucose Control With Metformin On 
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Complications In Patients With Type 2 Diabetes (UKPDS 34). Lancet, 352: 837-853, 
1998), have repeatedly demonstrated that the most effective way to prevent the long 
term complications of diabetes is by strictly maintaining blood glucose (BG) levels 
within a normal range using intensive insulin thet^y. 
S However, the same studies have also documented some adverse effects of 

intensive insulin therapy, the most acute of which is the increased risk of frequent 
severe hypoglycemia (SH), a condition defined as an episode of neuroglycopenia 
which precludes self-treatment and requires external help for recovery {See DCCT 
Research Group: Epidemiology of Severe Hypoglycemia In The Diabetes Control and 

10 Complications Trial. American Journal of Medicine^ 90: 4S0-4S9, 1991, and DCCT 
Research Group: Hypoglycemia in the Diabetes Control and Complications Trial. 
Diabetes^ 46: 271-286, 1997). Since SH can result in accidents, coma, and even death, 
patients and health care providers are discouraged from pursuing intensive therapy. 
Consequentiy, hypoglycraiia has been identified as a major bamer to improved 

15 glycemic control (Cryer PE: Hypoglycemia is the Limiting Factor in the Management 
Of Diabetes. Diabetes Metab Res Rev, 15: 42-46, 1999). 

Thus, patients with diabetes face a life-long optimization problem of 
maintaining strict glycemic control without increasing their risk of hypoglycemia. A 
major challenge related to this problem is the creation of simple and reliable methods 

20 that are capable of evaluating both patients' glycemic control and their risk of 
hypoglycemia, and that can be applied in their eveiyday environments. 

It has been well known for more than twenty years that glycosylated 
hemoglobin is a marker for tiie glycemic control of individuals with Diabetes Mellitus 
(Type I or Type II). Niunerous researches have investigated this relationship and 

25 have foimd that glycosylated heinoglobin generally reflects the average BG levels of a 
patient over the previous two months. Since in the majority of patients with diabetes 
the BG levels fluctuate considerably over time, it was suggested that the real 
connection between integrated glucose control and HbAu would be observed only in 
patients known to be in stable glucose control over a long period of time. 

30 Early studies of such patients produced an almost deterministic relationship 

between the average BG level in the preceding 5 weeks and HbAu, and this curvilinear 
association yielded a correlation coefGcient of 0.98 (See Aaby Svendsen P, Lauritzen 
T, Soegard U, Nerup J (1982). Glycosylated Hemoglobin and Steady-State Mean 
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Blood Glucose Concentration in Type 1 (Insulin-Dependent) Diabetes, Diabetoloria. 
23, 403-405). In 1993 the DCCT concluded that HbAic was the "logical nominee" for 
a gold-standard glycosylated hemoglobin assay, and the DCCT established a linear 
relationship between the preceding mean BG and UbAuiSee Santiago JV (1993). 
5 Lessons from the Diabetes Control and Complications Trial, Diabetes, 42, 1S49-1SS4). 
Guidelmes were developed indicating that an HbAie of 7% corresponds to a 
mean BG of 8.3 mM (ISO mg/dl), an HbAu of 9% corresponds to a mean BG of 1 1.7 
mM (210 mg/dl), and a 1% increase in HbAu corresponds to an increase in mean BG 
of 1 .7 mM (30 mg/dl, 2). The DCCT also suggested that because measuring the mean 
ID BG direcdy is not practical, one could assess a patient's glycemic control with a single, 
simple test, namely HbAu. However, studies clearly demonstrate that HbAu is not 
sensitive to hypoglycemia. 

Indeed, there is no reliable predictor of a patient's inunediate risk of SH from 
any data. The DCCT concluded that only about 8% of future SH could be predicted 
1 5 from known variables such as the history of SH, low HbAtc. and hypoglycemia 

unawareness. One recent review details the current clinical status of this problem, and 
provides options for preventing SH, that are available to patients and their health care 
providers (See BoUi, GB: How To Ameliorate The Problem of Hypoglycemia In 
Intensive As Well As Nonintensive Treatment Of Type I Diabetes. Diabetes Care^ 22, 
20 Supplement 2: B43-BS2, 1999). 

Contemporary home BG monitors provide the means for fi^uent BG 
measurements through Self-Monitoring of BG (SMBG). However, the problem with 
SMBG is that there is a missing link between the data collected by the BG monitors, 
and HbAic and hypoglycemia. In other words, there are currently no reliable methods 
25 for evaluating HbAu and recognizing imminent hypoglycemia based on SMBG 
readings (See Bremer T and Gough DA: Is blood glucose predictable from previous 
values? A solicitation for data. Diabetes 48:445-451, 1999). 

Thus, an object of this invention is to provide this missing link by proposing 
three distinct, but compatible, algorithms for evaluating HbAu and the risk of 
30 hypoglycemia from SMBG data, to be used to predict the short-term and long-temi 
risks of hypoglycemia, and the long-term risk of hyperglycemia. 

The inventors have previously reported that one reason for a missing link 
between the routinely available SMBG data and the evaluation of HbAu the .risk 
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of hypoglycemia, is that the sophisticated methods of data collection and clinical 
assessment used in diabetes research, are infirequently supported by diabetes-specific 
and mathematically sophisticated statistical procedures. 

Responding to the need for statistical analyses that take mto account the 
S specific distribution of BG data, the inventors developed a synmietrizing 

transformation of the blood glucose measurement scale {See Kovatchev BP, Cox DJ, 
Gonder-Frederick LA and WL Clarice (1997). Symmetization of the Blood Glucose 
Measurement Scale and Its Aiq)Iications, Diabetes Care. 20, 16SS-16S8) that worics as 
the follows. The BG levels are measured in mg/dl in the United States, and in mmol/L 
10 (or mM) in most other countries. The two scales are directly related by 18 mg/dl - 1 
mM. The entire BG range is given in most references as 1.1 to 33.3 mM, and this is 
considered to cover practically all observed values. According to the 
recommendations of tfie DCCT {See DCCT Research Group (1993) The Effect Of 
Intensive Treatment of Diabetes On the Development and Progression of Long-Term 
15 Complications of Insulin-Dependent Diabetes Mellitiis. New England Journal of 
Medicine^ 329, pp 978-986) the target BG range - also known as the euglycemic 
range- for a person witii diabetes is 3.9 to 10 mM, hypoglycemia occurs when the BG 
falls below 3.9 mM, and hyperglycemia is when the BG rises above 1 0 mM. 
Unfortunately, this scale is numerically asynunetric - the hyperglycemic range (10 to 
20 33.3mM) is wider tium tiie hypoglycemic range (1 . 1 to 3.9mM), and the euglycemic 
range (3.9 to lOmM) is not centered within the scale. The inventors correct this 
asymmetry by mtroducing a transformation,/(BG;, vAiich is a continuous fimction 
defined on die BG range [L 1, 33.3]^ having the two-parameter analytical form: 
f(BG.a,^^[(ln(BG)f-fil a,fi>0 
23 and which satisfies the assumptions: 

Al: f (33,1a, fi) = -/(Ha, fi) and 
A2: fOOAa^fi) = ./r5.P,a; py 
Next, fC) is multiplied by a third scaling parameter to fix the minimum and 
maximum values of the transformed BG range at - VlO and VIo respectively. These 
30 values are convenient since a random variable with a standard normal distribution has 
99.8% of its values witiiin tiie interval [- VlO , VlO ] . If BG is measured in mmol/1, 
when solved numerically with respect to the assumptions Al and A2, the parameters of 



4 



wo 01/72208 PCTAJSOl/09884 
the fimction/^BG, a,p)^a^ L026, fi = L861 and the scaUng parameter is y = 
1.794. IfBG is measured in mg/dl instead, the parameters are computed to be a = 
1.084. fi^ 5.381, and y = 1.509. 

Thus, when BG is measured in nunol/l, the synmietrizing transformation is 
5 f(BG) = 7. fSG;/-^^. 7.«tf77. and when BG is measured in mg/dl the 
symmetrizmg transformation \sf(BG) = 1.509[(ln (BG)/'^ -5.381]. 

On the basis of the symmetrizing transformation/(C^ the inventors introduced 
the Low BG Index - a new measure for assessing the risk of hypoglycemia &om 
SMBO readings (See Cox DJ, Kovatchev BP, Julian DM, Gonder-Frederick LA, 
10 Polonsky WH, Schlundt DG, Clarke WL: Frequency of Severe Hypoglycemia In 
IDDM Can Be Predicted From Self-Monitoring Blood Glucose Data. Journal of 
Clinical Endocrinology and Metabolism, 79: 1659-1662, 1994, and Kovatchev BP, 
Cox DJ, Gonder-FredCTick LA Young-Hyman D, Schlundt D, Clarke WL. Assessment 
of Risk for Severe Hypoglycemia Among Adults With IDDM: Validation of the Low 
15 Blood Glucose Index, Diabetes Care 21:1 870-1 875, 1998). Given a series of SMBG 
data the Low BG Index is computed as the average of lO./fBG/ taken for values of 
f(BG) <0 and 0 otherwise. Also suggested was a High BG Index, computed in a 
symmetrical to the Low BG Index manner, however this index did not find its practical 
application. 

20 Using the Low BG Index in a regression model the inventors were able to 

account for 40% of the variance of SH episodes in the subsequent 6 months based on 
the SH history and SMBG data, and later to enhance this prediction to 46% {See 
Kovatchev BP, Straume M, Farhi LS, Cox DJ: Estimating the Speed of Blood Glucose 
Transitions and its Relationship With Severe Hypoglycemia. Diabetes, 48: 

25 Supplement 1 , A363, 1 999). 

In addition, the inventors developed some data regarding HbAu and SMBG 
{See Kovatchev BP, Cox DJ, Straume M, Farhy LS. Association of Self-monitoring 
Blood Glucose Profiles with Glycosylated Hemoglobin. In: Methods in Enzvmoloev, 
vol. 321 : Numerical Computer Methods, Part C, Michael Johnson and Ludvig Brand, 

30 Eds., Academic Press, NY; 2000). 

These developments became a part of the theoretical background of this 
invention. In order to bring this theory into practice, several key theoretical 
components, among other things, as described in the following sections, were added. 
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In particular, three methods were developed for employing the evaluation of HbAic, 
long-term and short-term risk for hypoglycemia The development of these methods 
was, but not limited thereto, based on detailed analysis of data for 867 individuals with 
diabetes that included more than 300,000 SMBG readings, records of severe 
5 hypoglycemia and determinations of HbAu. 

The inventors have therefore sought to improve upon the aforementioned 
limitations associated with the conventional methods, and thereby provide simple and 
reliable methods that are capable of evaluating both patients' glycemic control and 
their risk of hypoglycemia, and that can be applied in their everyday environments. 

10 

SUMMARY OF THE INVENTIQN 

The invention includes a data analysis method and computer-based system for 
the simultaneous evaluation, from routinely collected SMBG data, of the two most 
important components of glycemic control in diabetes: HbAu and the risk of 

15 hypoglycemia. For the purposes of this document, self-monitoring of BG (SMBG) is 
defined as any method for determination of blood glucose at diabetic patients* natural 
envirorunent and includes the methods used by contemporary SMBG devices 
customarily storing 200-250 BG readings, as well as methods used by emerging 
continuous monitoring technologies. Given this broad definition of SMBG, this 

20 invention pertains directly to the enhancement of existing home blood glucose 

monitoring devices by introducing an intelligent data interpretation component capable 
of predictmg both HbAu and periods of increased risk of hypoglycemia, as well as to 
enhancement of future continuous monitoring devices by the same features. 

One aspect of the invention includes a method, system, and computer program 

25 product for evaluating HbAu from a predetermined period of collected SMBG data, 
for example 4-6 weeks. In one embodiment, the invention provides a computerized 
method and system for evaluating the HbAu of a patient based on BG data collected 
over a predetermined duration. The method includes computing weighted deviation 
toward high blood glucose (WR) and estimated rate of change of blood glucose (Dr) 

* 30 based on the collected BG data; estimating HbAu using a predetermined mathematical 
formula based on the computed WR and Dr, and providing a predetermined confidence 
interval for classification of said estimated value of HbAu* 
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Another aspect of the invention includes a method, system, and computer 
program product for estimating the long-term probability for severe hypoglycemia 
(SH). This method uses SMBG readings fiom a predetermined period, for example 4- 
6 weeks, and predicts the risk of SH within the following 6 months. In one 
5 embodiment, the invention provides a computerized method and system for evaluating 
the long term probability for severe hypoglycemia (SH) of a patient based on BG data 
collected over a predetermined duration* The method includes: computing weighted 
deviation toward low blood glucose (WL) and estimated rate of fall of blood glucose in 
the low BG range (DrDn) based on the collected BG data; estimating the number of 

10 future SH episodes using a predetermined mathematical formula based on the 

computed WL and DrDn; and defining a probability of incurring a select number of 
SH episodes respective to said estimated SH episodes. 

Still yet another aspect of the invention includes a method, system, and 
computer program product for identifying 24-hour periods (or other select periods) of 

15 increased risk of hypoglycemia. This is accomplished through the computation of the 
short-tenn risk of hypoglycemia using SMBG readings collected over the previous 24 
hours. In one embodiment, the invention provides a computerized method and system 
for evaluating the short term risk for severe hypoglycemia (SH) of a patient based on 
BG data collected over a predetermined duration. The method includes: computing 

20 weighted deviation toward low blood glucose (WL); determining Max(wl) by 
calculating maximum value of wl(BG;2)\ detemiining risk value by taking the 
geometric mean of WL and Max(wl) over the predetermmed duration; providing a 
predetermined threshold risk value; and comparing the determined risk value to the 
threshold risk value. 

25 These three aspects of the invention can be integrated together to provide 

continuous information about the glycemic control of an individual with diabetes, and 
enhanced monitoring of the risk of hypoglycemia. 

These and other objects, along with advantages and features of the invention 
disclosed herein, will be made more c^parent from the description, drawings and 

30 claims that follow. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other objects, features and advantages of the present 
invention, as well as the invention itself, will be more fiilly understood from the 
following description of preferred embodiments, when read togedier with the 
5 accompanying drawings in which: 

FIG. 1 is a flow chart illustrating the method of calculating the estimated HbAu 
and predicted HbAu confidence intervals in accordance with the present invention. 

FIG. 2 is a flow chart illustrating the method of calculating the estimated 
number of future SH episodes and the associated probability thereof in accordance 
10 with the present invention. 

FIG. 3 is a flow chart illustrating the method of calculating the estimated short 
term risk of an incurring imminent SH in accordance with the present invention. 

FIG. 4 is gr^>hical representation of a typical BG disturbances observed before 
and after an episode of severe hypoglycemia. 
15 FIG. S illustrates the action of the method for predicting short-term SH by 

presenting 10 weeks of data for Subject A (upper panel) and Subject B (lower panel). 
SH episodes are marked by triangle; a black line presents the risk value. When the risk 
threshold is crossed, the method indicates a subsequent high-risk period (gray bar). 

FIG. 6 is a functional block diagram for a computer system for implementation 
20 of the present invention. 

FIGS. 7-9 are schematic block diagrams of alternative variations of the present 
invention related processors, communication links, and isystems. 

DETAILED DESCRIPTION OF THE INVENTION 

25 The invention makes possible, but not lunited thereto, the creation of precise 

methods for the evaluation of diabetics* glycemic control, and include, firmware and 
software code to be used in computing the key components of the method. The 
inventive methods for evaluating HbAu, the long-term probability of SH, and the 
short-term risk of hypoglycemia, are also validated based on the extensive data 

30 collected, as will be discussed later in this document. Finally, the aspects of these 
methods can be combined in structured display or matrix. 
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Statiftnary pleasures of BG deviation 
According to the inventors' theory of BG symmetrization (See Kovatchev BP, 
Straume M, Cox DJ, Farhi LS. Risk Analysis of Blood Glucose Data: A Quantitative 
Approach to Optimizing the Control of Insulin Dependent Diabetes. J of Theoretical 
s Medicine, 3:1-10^2001) the natural clinical center of the BG measurement scale is at 
BG level of 1 1 2.S mg/dl (6.2S mmol/1) - a safe euglycemic value for a diabetes 
patient. 

Given this clinical center of the BG scale, the weighted deviations to the left 
(towards hypoglycemia) or to the right (towards hyperglycemia) are computed. The 
10 degree of weighting of these deviations will be rq)resented by parameters a and 6 
respectively as follows: 

wl(BG;a) = 10/(BGf if /(BG;<0 and 0 otherwise, and 
wr(BG;b) = 10.f(BGf if f(BG)>0 and 0 otherwise, 
where/(SG; is the BG symmetrization function presented in the backgroimd section. 
IS The weighting parameter a and 6 could be different, or the same for the left and right 
deviations. The inventors' data analyses demonstrated that the optimal for practical 
application parameto: values are a^2 (\^ch is the parameter value used for 
computation of the Low BG Index) and h^L Given a series of BG readings x/, ... 
x„, the average weighted deviations to the left and to the right of the clinical center of 
20 the BG scale are defined as: 

WL = -Y.^Kx,\2) and ira = -Y HT(x,a) respectively. 

These two measures of BG deviation do not depend on the timing of the BG readings, 
and therefore are stationary. In order to cloture the dynamics of BG change, measures 
of the BG rate of change are introduced as provided below. 

25 

Computation of BG risk rate nf chanpe 
Let xu X2, ... X;, be n SMBG readings of a subject recorded at time points 
tu t2, ... /«. This data is next transformed by calculating the numbers f(xi),f(x^) 
f( Xn) and draw a cubic spline S(t) passing through the points (tij(xi)), (t2f(x2,)) 
30 ( tnj( Xn)). Thus, thc function S(t) is a continuous function defined on the whole 
interval [ti , tj and such that S(tj)=f(xj\ fory=/, ...,it Also calculated are the set of 
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numbers Sk ^10.S(k'^ tif for t=OJ, ..., t„ - //. thus getting interpolated values at one- 
hour increments. 

Next, consider all couples of numbers Sii with consecutive indices: 

Co=(S(hSi), Q=(si,s^ , C2=(szS3)y and denote by Aftg,ibR set of all couples Ct, 

5 such that Sk > 5*+/ and by A/a>i the set of all couples C*, such that Sk < 5*+/ . 

Finally, let DrDn be the average of the numbers - st, provided that Q 6 
A/ah. and ZV be the average of the numbers j^+y - j*. provided that Q e A/,p+ Mdn- 

The numbers DrDn and Dr provide a measure for the rate of change of BG in a 
'*risk space " e.g. the fate of change of the risk associated with any BG level change. 
10 In addition, DrDn measures the rate of BG change only when BG goes down, i.e. 
DrDn evaluates how quickly the risk could increase when BG falls, while Z>r is a 
measiue of the overall vukierability of BG to fluctuations. It is further asserted that 
DrDn will be associated with risk for hypoglycemia (if someone's blood glucose could 
fall quickly, his/her risk for hypoglycemia would be higher), while Dr will be 
15 associated with the overall stability of BG. 



Software code (presented in SPSS control language^ 
The first is for vAien the BG readings are in mmol/L, and in this case the 
variable is BGMM. The second is for v^en the BG readmgs are in mg/dl, and m this 
20 case the variable is BGMG. 

If BG is measured in mmol/L, each BG reading is first transformed as follows: 

SCALEl=(ln{BGMM))**1.026 - 1.861 
RISK1=32 . 185*SCALE1*SCALE1 

If BG is measured in mg/dl, each BG reading is first transformed as follows: 

SCALE2- (In (BGMG) )**1. 08405 - 5.381 
RISK2-22 . 765*SCALE2*SCALE2 

Further, the left and right weighted deviations are computed as follows: 

WL=OWL=0 

IF (SCALE! le 0.0) WL»RISK1 
30 WR»0 

IF (SCALEl gt 0.0) WR=sqrt (RISKl) 
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Provided that the BG readings are equally spaced in time, or are interpolated at 
one-hour increments, the BG rate of change is computed as: 
Dr = RISK1(BG)-RISK1(BG-1) 
DrDn=0 

5 IF (SCALE le 0.0 and Dr gt 0) DrDn»Dr 

Finally, an aggregation pass through all BG readings for a subject will produce: 

WL = mean(WL) 
WR » mean{WH) 

Dr » mean(Dr), and DrDn » mean(DrDn) 

10 

Method for the Evaluation of HbA i. 
A preferred embodunent of HbAu evaluation method 100 according to the 
invention is illustrated in FIG. 1 . In a first step 1 02, SMBG data is collected over a 
predetermined period of time. For example, the SMBG data is collected ov^ 4-6 

IS weeks with a frequency of 3-S BG measurements per day, of which are transformed by 
the code or formulas presented in the previous section. DifiTerent formulas are to be 
used iftheBG measurements are stored in mg/dl, or in mmol/1. One skilled in the art 
would appreciate that various levels, durations, and frequencies can be employed, bi a 
step 104, weighted deviation towards high blood glucose (WR) and estimated rate of 

20 change of blood glucose (Dr) is computed using the formula / code discussed above. 
In a step 106, an estimate of HbAu from self-monitoring data is computed using the 
linear function: EstHBAlc = 0.9008«WR - 0.8207*DR + 6.7489. It is noted that the 
coefBcients of this function are derived from data for 867 individuals with diabetes, 
and one would recognize that further data accumulation may update these coefficients. 

25 In step 1 08 HbAu estimate categories representing a range of values for estimated 
HbAu are defined according to Table 1. 



Table I : Defining categories on the basis of EstHBAlc: 



EstHBAlc 


<7.8 


7.8-8.5 


8.5-9.0 


9.0-9.6 


9.6-10.3 


10.3-11.0 


> 11.0 


Category 


1 


2 


3 


4 


5 


6 


7 
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In Step 1 10 predicted confidence intervals for corresponding HbAu are derived 
according to Table 2. 



Table 2: Predicted 95% confidence intervals for classification of HbAu: 



Category 


1 


2 


3 


4 


5 


6 


7 


HBAie 


<8.0 


8.0-8.5 


8.5-9.0 


9.0-9.5 


9.5-10.1 


10.1-11.0 


>11.0 



5 In step 1 12, the estimated HbAu fiom step 106 is assigned in one of the 

categories provided in Table 1 and / or Table 2. 



Empirical Validation of Evaluat ion of HbA t. 
The intervals for HbAic in Table 2 are based on extensive research. To validate 

10 these intervals we analyzed SMBG and HbAic data fix)m 867 subjects with diabetes. 
All subjects were instructed to use BG memory meters for six months and to measure 
their BG two to four times a day. During the same period S to 8 HbAic assays were 
performed for each subject The memory meter data were electronically downloaded 
and stored in a computer for further analysis. This procedure produced a database 

15 containing more than 300,000 SMBG readings and 4,180 HbAu assays taken over six 
months. Analysis of variance was conducted to compare HbAu in the seven categories 
identified in Table 1. The five categories were highly significantly different, with 
F=91 and p<0.00001. Moreover, the average HbAu was significantly different for 
each pair of categories as demonstrated by Duncan's ranges, with pO.Ol . 

20 Also, 95% confidence intervals were computed for the mean value of HbAu in 

each of the seven categories. These confidence intervals were used as a basis for 
computing the HbAu intervals presented in Table 2. Post-hoc analysis of the 
classification power of this method demonstrated that the method was well protected 
against extreme errors such as incorrectly classifying HbAu in category 1 , 2 or 3 on 

23 the basis of SMBG while the actual HbAu was greater than 9.5%, or classifying HbAu 
in category 5, 6 or 7 vAdle ttie actual HbAu was below 9.0%. 

In sununary, after an initial 4-6 weeks.of SMBG readings the computerized 
method computes an interval estimate for the value of HbAu that can be used to track 
patients' changes in glycemic control in the high BG range. 
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Method for Evaluation of the Long-Term Probability for Severe Hypoglycemia (SH) 
A preferred embodiment of long-term probability for SH evaluation method 
200 according to the invention is illustrated in FIG. 2. In a first step 202, SMBG data 
is collected over a predetermined period of time. For example, the SMBG data is 
5 collected over 4*6 weeks with a frequency of 3-5 BO measurements per day, of which 
are transformed by the code or formulas presented immediately above. Different 
formulas are to be used if the BG measurements are stored in mg/dl, or in nunol/l; 
One skilled in the art would appreciate that various levels, durations, and fi^uencies 
can be employed In a step 204, WL and DrDn are computed using the formula / code 
10 as discussed above. In step 206, an estimate of the number of future SH episodes is 
computed using the linear function: 

EstNSH = 3.3613*WL - 4.3427*DrDn - 1.2716. 
It is noted that the coefBcients of this function are derived from data for 1 81 
individuals with diabetes, and one would t^preciate that further data accumulation 
1 5 may update these coefficients. It is further noted that this formula provides a single 
value estimate for the number of future SH episodes and that through additional 
methodologies, as discussed below, categories are provided with ranges and 
confidence levels for enhanced clinical applications. In step 208, estimated number of 
SH episodes (estNSH) categories representing a range of values for estNSH are 
20 defined according to Table 3 . 

Table 3: Classification of EstNSH: 



EstNSH 


< 0.775 


0.775-3.750 


3.750 - 7.000 


> 7.000 


Category 


1 


2 


3 


4 



25 



In step 210, respective to the estNSH categories, the probability of incurring 0, 
1-2, or more than 2 SH episodes in the following six months is derived, as represented 
in table 4. 





Category 1 


Category 2 


Category 3 


Category 4 


OSH 


90% 


50% 


25% 


<20% 


1-2 SH 


10% 


25% 


25% 


>2SH 


25% 


50% 


> 80% 
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In Step 212, the EstNSH from step 206 is assigned in one of the categories 
provided in Table 3 and / or Table 4. 



Empirical Validation of Evaluation of the Long-Term Probability for SH 
5 One-hundred-eighty-one adults with Type 1 diabetes (mean age 37 years, 

duration of diabetes 1 8 years) used memory meters to collect more than 34,000 SMBG 
over a month. Then for the next six months they recorded in diaries any occurrence of 
SH. The SMBG data were mathematically transformed and an a linear regression 
model was used to predict future severe hypoglycemia resulting in a highly significant 

10 model (F=36.3, p<0.0001) and multiple R of 55%. 

All subjects were classified into 4 categories using the present long-term SH 
method. The average number of future SH episodes in categories 1, 2, 3, and 4 was 
0.3, 2.0, S.O, and 9.75 respectively. Analysis of variance demonstrated highly 
significant diflferences between these categories, F=19.0, p<0.0001. 

15 In summary, a linear combination of the Low BG Index and the rate of drop of 

BG as measured in "risk space" provide an accurate assessment of the long-term risk 
of SH. Because it is based on SMBG records that are automatically stored by many 
reflectance meters, this is an effective and clinically useful indicator of patients' 
glycemic control in the low BG range. 

20 

Method for the Evaluation of the Sho rt-term fwithin 24 hours) Risk of Hvpoclvcemia 

A preferred embodiment of short term risk of SH evaluation method 300 
according to the invention is illustrated in FIG. 3. In a first step 302, SMBG is data is 
collected over a predetermined short term period. For example, the SMBG data is 

25 collected over a 24 hour period, with a fi«quency of 3-5 BG measurements per day - 4 
or more readmgs, as a nominal level according to data analyses. One skilled in the art 
would appreciate that various levels, periods (durations), and firequencies can be 
employed. In a step 304 WL(24) and Max(wl) is computed bom all readings collected 
within the preceding 24 hours, wherein the maxunum value of wl(BG;2) is Max(wl). 

30 bi step 306, the risk value is by taking the geometric mean of WL and Max(wl) over 
the 24 hour period, >^erein said risk value is mathematically defined as: 

Risk(24) = V^I(24)«A/ar(w/); 
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In Step 308 a threshold risk value is detennined. In step 3 1 0 the estimated risk 
value is compared to the threshold risk value. For example, if the threshold risk value 
is set at 1 7, then if Risk(24)>l 7, then-based on the SMBG data collected over the 
previous 24 hours- the resultant indication is a high risk of the patient incurring 
5 imminent hypoglycemia. In other words, this is a decision-making rule that considers 
a 24-hour period of SMBG data and judges whether this period is likely to precede an 
imminent hypoglycemia episode. The threshold value of 17 is derived fix>m an 
extensive data set, however, it is recognized that it is possible that this value maybe 
adjusted with further accumulation of data or for additional objectives. 

10 

Empirical Validation of Evaluation of the Short-term Risk of Hvpoglvcemia: 
Eighty-five individuals were recmited through advertisement in newsletters, 
diabetes clinics, and through direct referrals. The inclusion criteria were: 1) age of 21- 
60 years; 2) type I diabetes with at least two years duration, and insulin use since the 

1 5 time of diagnosis; 3) at least 2 documented SH episodes in the past yeai^ and 4) routine 
use of SMBG devices for diabetes monitoring. The participants were instructed to use 
the meter 3-S times a day, and to record in monthly diaries any SH episodes, including 
the exact dates and times of their occurrences. SH was defined as severe 
neuroglycopenia that results in stupor or unconsciousness and precludes self-treatment. 

20 For each subject the study continued 6-8 months and each month the subject's meter 
was downloaded and tiie SH diary was collected. The memory capacity of the meters 
was sufficient, and the downloading was often enough, so that no BG data were lost. 
No changes were made in the participants' diabetes management routine, nor were any 
additional treatments administered during the study. 

25 During the study a total of 7S,49S SMBG readings (on average 4.0±1 .5 per 

subject per day) were downloaded fiom the participants' memory meters, and 399 
(4.7±6.0 per subject) SH episodes were recorded in their diaries. An important 
finding, among other things, was that episodes of moderate or severe hypoglycemia are 
preceded and followed by measurable BG disturbances. In the 24-hour period before 

30 an SH episode the Low BG Index (e.g. WL) rose (p<0.001), the average BG was lower 
(p=0.001), and the BG variance increased (p=0.001). In the 24 hours following the SH 
episode, the Low BG Index and BG variance remained elevated (p<0.001), but the 
average BG returned to its baseline. 
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To this end» FIG. 4 is graphical representation of a typical BG disturbance 
observed before and after an episode of severe hypoglycemia. In the period 48 to 24 
hours before the SH episode, the average BG level decreased and the variance of BG 
increased. In the 24-hour period immediately preceding the SH episode, the average 
5 BG level dropped further and the variance of BG continued to increase. In die 24-hour 
period following the SH episode, the average BG level normalized, but the BG 
variance remained greatly increased. Both the average BG and its variance returned to 
their baseline levels within 48 hours after the SH episode. 

As such, as part of tiie invention, the disturbances presented in FIG. 4 are 

10 quantified from SMBG data to enable the evaluation of the short-tenn risk of 

hypoglycemia. The cutoff value of Risk(24)°l 7 is derived from an optunization along 
the following restrictions: 1) the method had to predict a maximum percentage of SH 
episodes, i.e. to identify as risky a maximum percentage of 24-hour periods pieceding 
SH, and 2) to prevent overestimadon of the risk, the method had to identify as risky no 

15 more that 15% of the total time of the study (one day a week on average). The cutoff 
risk value of 1 7 was held constant for all subjects. The reason for choosing the value 
of 1 S% was to prevent the patients from becoming irritated with an overabundance of 
'Yalse alarms" and then ignoring 'true alarms." In practice, a patient's physician can 
select an alternate value depending on the severity of the patient's diabetes and 

20 particular objectives. 

The following example illustrates the action of the algorithm on the SMBG 
data of two participants in the study. FIG. S presents ten weeks of data for Subject A 
(upper panel) and Subject B (lower panel). SH episodes are marked by triangles; a 
black curve presents the risk value. When the risk threshold (the horizontal line at 

25 Risk=l 7) is crossed, the algorithm indicates a subsequent high-risk period (gray bar). 
For Subject A, 7 out of 9 SH episodes are predicted and there are 5 false alarms, e.g. 
high-risk periods that did not result in SH; for Subject B there are 3 false alarms and 
the only SH episode is predicted. It is obvious that Subject B's risk values vAiea 
compared to Subject A*s risk values, include more and higher deviations. For both 

30 subjects, all SH episodes were accompanied by supercritical risk values, and about half 
of all large deviations were accompanied by one or more SH episode. 

Across all participants in the study, 44% of all recorded SH episodes were 
preceded, within 24 hours, by a high-risk period, and 50% were preceded, within 48 
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hours, by a high-risk period. If only periods with either at least 3, or at least 4 SMBG 
measurements were considered, the accuracy of the latter prediction increased to 53% 
and 57%, respectively. Post-hoc analysis of BG levels occurring during, or 
immediately after, high-risk periods that were not followed by an SH episode, i.e. 
5 during, or inmiediately after false alarms, demonstrated that the average per subject 
minimum of such BG levels was 2.3±0.2 nunol/l versus 5.9±1.7 mmol/1 (t=19.5, 
pO.OOOI) for all non-risk periods, including all SH episodes that remained 
unaccounted for. This indicates that, although symptomatic SH did not occur, BG 
levels following high-risk periods were notably low. 
10 In summary, the inventors simulated the action of the short-term risk method 

on a 6-month series of SMBG readings for 8S individuals with Type I diabetes. With 
four or more SMBG readings per day, at least 50% of all episodes of SH could be 
anticipated. Even when symptomatic SH did not occur, the algorithm predicted 
episodes of moderate hypoglycemia. 

15 

Integration of the Three Methods 
The three methods of this mvention, as discussed above and illustrated in FIGS. 
1-3, utilize the same series of SMBG data. Therefore, fix)m an SMBG-device point of 
view, a unified display or matrix of the results of these three methods could be made 
20 similar to the grid output presented below: 





EstHBA categories (Algorithm 1) 


1 


2 


3 


4 


5 


6 


7 


EstNSH 
categories 


1 
















2 
















3 
















4 

















25 



Thus, for example, the output for subject 1 (Ss 1) shown in the above grid 
indicates that this person is likely to have HbAu between 9 and 9.5%, and has a 90% 
chance not to experience severe hypoglycemia in the subsequent 6 months. The output 
for subject 2 (Ss 2) indicates that this person is likely to have HbAu below 8%, and 
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has a greater than 80% chance to experience at least 3 SH episodes in the subsequent 6 
months. 

In addition to this grid-ou^ut, the short term risk method provides a continuous 
tracking of the risk of imminent hypoglycemia and can be used to sound an alarm 
5 ivhen this risk becomes high. 

The method of the invention may be implemented using hardware, software or 
a combination thereof and may be implemented in one or more computer systems or 
other processing systems, such as personal digit assistants (PDAs). In an example 
embodiment, the invention was implemented in software running on a general purpose 

10 computer 900 as illustrated in FIG. 6. Computer system 600 includes one or more 
processors, such as processor 604. Processor 604 is connected to a communication 
infrastructure 606 (e.g., a communications bus, cross-over bar, or network). Computer 
system 600 includes a display inter&ce 602 that forwards gnqphics, text, and other data 
ftom the communication infrastructure 606 (or from a frame buffer not shown) for 

1 5 display on the display unit 630. 

Computer system 600 also includes a main memory 608, preferably random 
access memory (RAM), and may also include a secondary memory 610. The 
secondary memory 610 may include, for example, a hard disk drive 612 and/or a 
removable storage drive 614, representing a floppy disk drive, a magnetic tape drive, 

20 an optical disk drive, etc. The removable storage drive 614 reads from and/or writes to 
a removable storage unit 618 in a well known manner. Removable storage unit 618, 
represents a floppy disk, magnetic tape, optical disk, etc. which is read by and written 
to by removable storage drive 614. As will be iq;>preciated, the removable storage unit 
6 1 8 includes a computer usable storage medium having stored therein computer 

25 software and/or data. 

In alternative embodiments, secondary memory 610 may include other means 
for allowing computer programs or other instructions to be loaded into computer 
system 600. Such means may include, for example, a removable storage xmit 622 and 
an interface 620. Examples of such removable storage units/interfrtces include a 

30 program cartridge and cartridge interface (such as that found in video game devices), a 
removable memory chip (such as a ROM, PROM, EPROM or EEPROM) and 
associated socket, and other removable storage units 622 and interfaces 620 which * 



18 



wo 01/72208 PCT/USOl/09884 
allow software and data to be transferred from the removable storage unit 622 to 
computer system 600. 

Computer system 600 may also include a communications interface 624. 
Communications interface 624 allows software and data to be transferred between 
S computer system 600 and external devices. Examples of commimications interface 624 
may include a modem, a network interface (such as an Ethernet card^ a 
conmiunications port, a PCMCIA slot and card, etc. Software and data transferred via 
communications inter&ce 624 are in the form of signals 628 v^ch may be electronic, 
electromagnetic, optical or other signals capable of being received by communications 
10 interface 624. Signals 628 are provided to communications interface 624 via a 

communications path (i.e., channel) 626. Channel 626 carries signals 628 and may be 
implemented using wire or cable, fiber optics, a phone line, a cellular phone link, an 
RF link and other conunimications channels. 

In this document, the terms "computer program mediimi" and "computer usable 
] s medium" are used to generally refer to media such as removable storage drive 914, a 
hard disk installed in hard disk drive 612, and signals 628. These computer program 
products are means for providing software to computer system 600. The invention 
includes such computer program products. 

Computer programs (also called computer control logic) are stored in main 
20 memory 608 and/or secondary memory 610. Computer programs may also be received 
via communications interfiice 624. Such computer programs, when executed, enable 
computer system 600 to perform the features of the present invention as discussed 
herein. In particular, the computer programs, when executed, enable processor 604 to 
perform the functions of the present invention. Accordingly, such computer programs 
25 represent controllers of computer system 600. 

In an embodiment vdiere the invention is implemented using software, the 
software may be stored in a computer program product and loaded into computer 
system 600 using removable storage drive 614, hard drive 612 or conununications 
interface 624. The control logic (software), vAica executed by the processor 604, 
30 causes the processor 604 to perform the functions of the invention as described herein. 
In another embodiment, the invention is implemented primarily in hardware 
using, for example, hardware components such as ^plication specific integrated 
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circuits (ASICs). Implementation of the hardware state machine to perform the 
functions described herein will be apparent to persons skilled in the relevant art(s). 

In yet another embodiment, the invention is implemented using a combination 
of both hardware and software. 
5 In an example software embodiment of the invention, the methods described 

above were implemented in SPSS control language, but could be implemented in other 
programs such as, but not limited to, C + + programming language. 

FIGS. 7-9 show block diagrammatic representation of alternative 
embodiments of the invention. Referring FIG. 7, there is shown a block diagrammatic 

10 representation of the system 710 essentially comprises the glucose meter 728 used by a 
patient 712 for recording, inter alia, insulin dosage readings and measured blood 
glucose C'BG") levels. Data obtamed by the glucose meter 728 is preferably 
transferred through q>propriate communication links 714 or data modem 732 to a 
processmg station or chip, such as a personal computer 740, PDA, or cellular 

1 5 telephone. For instance, data stored may be stored within the glucose meter 728 and 
may be directly downloaded into the personal computer 740 through an appropriate 
interface cable. An example is the ONE TOUCH monitoring system or meter by 
LifeScan, Inc. which is compatible with IN TOUCH software which includes an 
interface cable to down load the data to a personal computer. 

20 The glucose meter is common in the industry and includes essentially any 

device that can functions as a BG acquisition mechanism. The BG meter or 
acquisition mechanism, device, tool, or system includes various conventional methods 
directed toward drawing a blood sample (e.g. by fingerprick) for each test, and a 
determination of the glucose level using an instrument that reads glucose 

25 concentrations by electromechanical or claorimetric methods. Recently, various 
methods for determining tiie concentration of blood analytes without drawing blood 
have been developed. For example, U.S. Pat. No. 5,267,1 52 to Yang et al. describes a 
noninvasive technique of measuring blood glucose concentration using near-IR 
radiation diffuse-reflection laser spectroscopy. Similar near-IR spectrometric devices 

30 are also described in U.S. Pat. No. 5,086,229 to Rosenthal et al. and U.S. Pat No. 
4,975,581 to Robinson et al. 

U.S. Pat. No. 5,139,023 to Stanley describes a transdermal blood glucose 
monitoring apparatus that relies on a permeability enhancer (e.g., a bile salt) to 
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Militate transdermal movement of glucose along a concentration gradient established 
between interstitial fluid and a receiving medium. U.S. Pat. No. 5,036,861 to 
Sembrowich describes a passive glucose monitor that collects perspiration through a 
skin patch, where a cholinergic agent is used to stimulate perspiration secretion from 
S the eccrine sweat gland. Similar perspiration collection devices are described in U.S. 
Pat. No. 5,076,273 to Schoendorfer and U.S. Pat, No. 5,140,985 to Schroeder. 

In addition, U.S. Pat No. 5;279,543 to Glikfcld describes the use of 
iontophoresis to noninvasively sample a substance through skin into a receptacle on 
the skin surface. Glikfeld teaches that this sampling procedure can be coupled with a 

10 glucose-specific biosensor or glucose-specific electrodes in order to monitor blood 
glucose. Moreover, International Publication No. WO 96/001 10 to Tamada describes 
an iontophoretic apparatus for transdermal monitoring of a target substance, vdierein 
an iontophoretic electrode is used to move an analyte into a collection reservoir and a 
biosensor is used to detect the target analyte present in the reservoir. Finally, U.S. Pat. 

1 5 No. 6,144,869 to Bemer describes a sampling system for measuring the concentration 
of an analyte present.. 

Further yet, the BG meter or acquisition mechanism may include indwelling 
catheters and subcutaneous tissue fluid sampling. 

The computer or PDA 740 includes the software and hardware necessary to 

20 process, analyze and interpret the self-recorded diabetes patient data in accordance 
with predefined flow sequences (as described above in detail) and generate an 
appropriate data interpretation ou^ut Preferably, die results of the data analysis and 
interpretation performed upon the stored patient data by the computer 740 are 
displayed in the form of a paper report generated through a printer associated v^th the 

25 personal computer 740. Alternatively, the results of the data interpretation procedure 
may be directly displayed on a video display unit associated with the computer 740. 

FIG. 8 shows a block diagranmiatic representation of an alternative 
embodiment having a diabetes management system that is a patient-operated apparatus 
810 having a housing preferably sufficiently compact to enable apparatus 810 to be 

30 hand-held and carried by a patient A strip guide for receiving a blood glucose test strip 
(not shown) is located on a surface of housing 816. Test strip is for receiving a blood 
sample fiom the patient 812. The apparatus includes a microprocessor 822 and a 
memory 824 connected to microprocessor 822. Nficroprocessor 22 is designed to 
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execute a computer program stored in memory 824 to perform the various calculations 
and control functions as discussed in great detail above. A keypad 8 1 6 is connected to 
microprocessor 822 through a standard keypad decoder 826. Display 814 is connected 
to microprocessor 822 through a display driver 830. Microprocessor 822 
5 communicates with display driver 830 via an interface, and display driver 830 updates 
and refreshes display 814 under the control of microprocessor 822. Speaker 854 and a 
clock 856 are also connected to microprocessor 822. Speaker 854 operates under the 
control of microprocessor 822 to emit audible tones alerdng the patient to possible 
future hypoglycemia. Clock 856 supplies the current date and time to microprocessor 
10 822. 

Memory 824 also stores blood glucose values of the patient 812, the insulin 
dose values, the insulin types, and the parameter values used by microprocessor 822 to 
calculate future blood glucose values, supplemental insulin doses, and carbohydrate 
supplements. Each blood glucose value and insulin dose value is stored in memory 

1 S 824 with a corresponding date and time. Memory 824 is preferably a non-volatile 
memory, such as an electrically erasable read only memory (EEPROM). 

Apparatus 810 also includes a blood glucose meter 828 connected to 
microprocessor 822. Glucose meter 828 is designed to measure blood samples received 
on blood glucose test strips and to produce blood glucose values fiom measurements of 

20 the blood samples. As mentioned previously, such glucose meters are well known in 
the art. Glucose meter 828 is preferably of the type Mfidch produces digital values 
which are output directly to microprocessor 822. Alternatively, blood glucose meter 
828 may be of the type which produces analog values. In this alternative embodiment, 
blood glucose meter 828 is cormected to microprocessor 822 through an analog to 

25 digital converter (not shown). 

Apparatus 810 further includes an input/output port 834, preferably a serial 
port, which is connected to microprocessor 822. Port 834 is cormected to a modem 832 
by an mterface, preferably a standard RS232 inter&ce. Modem 832 is for establishing 
a conmiunication link between apparatus 810 and a personal computer 840 or a 

30 healthcare provider computer 838 through a communication networic 836. Specific 
techniques for connecting electronic devices through connection cords are well known 
in the art. Another alternative example is "bluetooth** technology communication. 
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Alternatively, FIG. 9 shows a block diagrammatic representation of an 
alternative embodiment having a diabetes management system that is a patient- 
operated apparatus 910, similar as shown in FIG. 8, having a housing preferably 
sufiGciently compact to enable the apparatus 910 to be hand-held and carried by a 
S patient. However, the present embodiment includes a separate or detachable glucose 
meter or BG acquisition mechanism 928. 

Accordingly, the embodiments described herein are capable of being 
implemented over data conmumication networks such as the internet, making 
evaluations, estimates, and information accessible to any processor or computer at any 
10 remote location, as depicted in FIGS. 6-9 and/or U.S. Pat No. 5,851,186 to Wood, of 
which is hereby incorporated by reference herein. Alternatively, patients located at 
remote locations may have the BG data transmitted to a central healthcare provider or 
residence, or a different remote location. 

In sununary, the invention proposes a data analysis computerized method and 
15 system for the simultaneous evaluation of the two most important components of 
glycemic control in individuals with diabetes: HbAu and the risk of hypoglycemia. 
The method, M^le using only routine SMBG data, provides, among other things, three 
sets of output 

The potential implementations of the method, system, and computer program 
20 product of the invention is that it provides the following advantages, but are not limited 
thereto. First, the invention enhances existing home BO monitoring devices by 
producing and displaying: 1) estimated categories for HbAu, 2) estimated probability 
for SH in the subsequent six months^ and 3) estimated short-term risk of hypoglycemia 
(i.e. for the next 24 hours). The latter may include warnings, such as an alarm, that 
25 indicates inuninent hypoglycemic episodes. These three components can also be 
integrated to provide continuous information about the glycemic control of individuals 
with diabetes, and to enhance the monitoring of their risk of hypoglycemia. 

As a second advantage, the invention enhances existing software or hardware 
that retrieves SMBG data. Such software or hardware is produced by virtually every 
30 manufecturer of home BG monitoring devices and is customarily used by patients and 
health care providers to interpret SMBG data. The methods and system of the 
invention can be directly incorporated into existing home blood glucose monitors, or 
used for the enhancement of software that retrieves SMBG data, by introducing a data 
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interpretation component capable of predicting both HbAic and periods of increased 
risk of hypoglycemia. 

Still yet another advantage, the invention evaluates the accuracy of home BG 
monitoring devices, both in the low and high BG ranges, and over the entke BG scale. 
5 Moreover, another advantage, the invention evaluates the effectiveness of 

various treatments for diabetes. 

Further still, as patients with diabetes &ce a life-long optimization problem of 
maintaining strict glycemic control without increasing their risk of hypoglycemia, die 
present invention alleviates this related problem by use of its simple and reliable 
10 methods, i.e., the invention is oqpable of evaluating both patients' glycemic control 
and their risk of hypoglycemia, and at the same time iq>plying it in their everyday 
environments. 

Additioiudly, the invention provides the missing link by proposing three 
distmct, but compatible, algorithms for evaluating HbAu and the risk of hypoglycemia 

15 fiom SMBG data, to be used to predict the short-term and long-term risks of 
hypoglycemia, and the long-term risk of hyperglycemia. 

Finally, another advantage, the invention evaluates the effectiveness of new 
insulin or insulin delivery devices. Any manufacturer or researcher of insulin or 
insulin delivery devices can utilize the embodiments of the mvention to test the relative 

20 success of proposed or tested insulin types or device delivery designs. 

The invention may be embodied in other specific forms without departing fix>m 
the spirit or essential characteristics thereof. The foregoing embodiments are therefore 
to be considered in ail respects illustrative rather than limiting of the invention 
described herein. Scope of the invention is thus indicated by the appended claims 

25 rather than by the foregomg description, and all changes which come within the 
meaning and range of equivalency of the claims are therefore intended to be embraced 
thereirt 
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CLAIMS 
We claim: 

1 1 . A computerized method for evaluating the HbAu of a patient based on BO 

2 data collected over a predetermined duration, said method comprising: 

3 computing weighted deviation toward high blood glucose (WR) and estimated rate of 

4 change of blood glucose (Dr) based on said collected BG data; and 

5 estimating HbAu using a predetermined mathematical formula based on said 

6 computed WR and Dr. 

1 2. The method of claim 1 , wherein: 

2 said computed WR is mathematically defined fi^om a series of BG readings x/, X2, ... Xn 

3 taken at time points //, r^, r„ as: 

W M 

5 ^ere: 

6 wr(BG;b) = 10,f(BG)^ if f(BG)>0 and 0 otherwise, 

7 b=l, representing a weighting parameter, and 
8 

9 said computed DR is mathematicaUy defined as: 

10 /)r = average of - 5*, 

11 where: 

12 5*=iO.Sf*+r/ forJt=0.i,...,/„.r/, 

13 S(tjH(xj\ for j^l n. 

1 3. The method of claim 1, wherein said the estimate of HbA]cfix>m said BG 

2 monitoring data is mathematically defined as : 

3 Estimated HbAic= 0.9008(WR) - 0.8207(DR) + 6.7489. 

1 4. The method of claim 1, fiirther comprising: 

2 defining predetermined categories for the estimate of HbAu , each of said HbAu 

3 estimate categories representing a range of values for estimated HbAic; and 

assigning said estimated HbAu to at least one of said HbAic estimate categories. 
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1 S. The method of claim 4, wherein said HbAic estimate categories are defined as 

2 follows: 

3 classified category 1 , herein said estimated HbAic is less than about 7.8; 

4 classified category 2, wherein said estimated HbAu is between about 7.8 and about 

5 8.5; 

6 classified category 3, i^erein said estimated HbAu is between about 8.5 and about 

7 9.0; 

8 classified category 4, M^erein said estimated HbAu is between about 9.0 and about 

9 9.6; 

10 classified category 5, \Aerein said estimated HbAu is between about 9.6 and about 

11 10.3; 

12 classified category 6. wherein said estimated HbAu is between about 1 0.3 and about 

13 11.0; and 

14 classified category 7, wherein said estimated HbAu is above about 1 1.0. 

1 6. The method of claim S, further comprising: 

2 definii^ predicted confidence intervals for corresponding said HbAu estimate 

3 categories, wherein said predicted confidence intervals are defined as follows: 

4 said classified category 1 corresponds with a predicted HbAu less than about 

5 8.0; 

6 said classified category 2 corresponds with a predicted HbAu between about 

7 8.0 and about 8.5; 

8 said classified category 3 correqx)nds with a predicted HbAu between about 

9 8.5 and about 9.0; 

10 said classified category 4 corresponds with a predicted HbAu between about 

11 9.0 and about 9.5; 

12 said classified category 5 corresponds with a predicted HbAu between about 

13 9.5 and about 10.1; 

14 said classified category 6 corresponds with a predicted HbAu between about 

15 10.1 andabout 11.0; and 
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said classified category 7 corresponds with a predicted HbAic above about 

1 7. The method of claim 4, further comprising: 

2 defining predicted confidence intervals for corresponding said HbAu , each of said 

3 predicted confidence intervals representing a range of values for HbAu. 

1 8. The method of claim 7» wherein said predicted HbAu confidence intervals 

2 have about a 95% confidence level. 



1 9. A computerized method for evaluating the HbAu of a patient based on BG 

2 data collected over a predetermined duration, said method comprising: 

3 computing weighted deviation toward high blood glucose (WR) and estimated rate of 

4 change of blood glucose (Dr) based on said collected BG data; 

5 estimating HbAu using a predetermined mathematical formula based on said 

6 computed WR and Dr; and 

7 providing a predetermined confidence interval for classification of said estimated 

8 value of HbAu. 

1 10. The method of claim 9, v^erein: 

2 said confidence interval is between about 83% to about 95%. 

1 1 1 . A system for evaluating HbAu of a patient based on BG data collected over a 

2 predetermined duration, said system comprising: 

3 a database component operative to maintflin a database identifying said BG data; 

4 a processor progranmied to: 

5 compute weighted deviation toward high blood glucose (WR) and estimated 

6 rate of change of blood glucose (Dr) based on said collected BG data; and 

7 estimate HbAu using a predetermined mathematical formula based on said 

8 computed WR and Dr. 
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1 12. The system of claim 1 1 , herein: 

2 said computed WR is mathematically defined bom a series of BG readings Xi, X2» ... 

3 Xn taken at time points th th tn as: 

4 WR^^Y.wr{x,^) 

5 where: 

6 wr(BG;b) - lO.ffBG/ if f(BG)>0 and 0 otherwise, 

7 b =1 , representing a weighting parameter, and 
8 

9 said computed OR is mathematically defined as: 

10 Dr = average of sm - 

1 1 where: 

12 Sk^lOSO^-^ tsf for k=0,l...,tn-ih 

13 S(tj)^f(xjl for J=l n. 

1 13. The system of claim 1 1, herein said the estimate of HbAie from said BG 

2 monitoring data is mathematically defined as : 

3 Estimated HbAu= 0.9008(WR) - a8207(DR) + 6,7489. 

1 14. The system of claim 1 1, wherein said processor being further programmed to: 

2 define predetermined categories for the estimate of HbAic , each of said HbAfc 

3 estimate categories representing a range of values for estimated HbAic; and 

4 assign said estimated HbAu to at least one of said HbAu estimate categories. 

1 IS. The system of claim 14, wherein said HbAic estimate categories are defined as 

2 follows: 

3 classified category 1 , wherein said estimated HbAu is less than about 7.8; 

4 classified category 2, \\iierein said estimated HbAu is between about 7.8 and about 

5 8.5; 

6 classified category 3, wherein said estimated HbAu is between about 8.S and about 

7 9.0; 
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8 classified category 4, wherein said estimated HbAu is between about 9.0 and about 

9 9.6; 

10 classified category S, wherein said estimated HbAu is between about 9.6 and about 

11 10.3; 

12 classified category 6, wherein said estimated HbAic is between about 1 0.3 and about 

13 11.0; and 

14 classified category 7, wherein said estimated HbAu is above about 1 1 .0. 

1 1 6. The system of claim 1 S, wherein said processor being further programmed to: 

2 define predicted confidence intervals for corresponding said HbAic estimate 

3 categories, wherein said predicted confidence intervals are defined as follows: 

4 said classified category 1 corresponds with a predicted HbAu less than about 

5 8.0; 

6 said classified category 2 corresponds with a predicted HbAu between about 

7 8.0 and about 8.S; 

8 said classified category 3 corresponds with a predicted HbAu between about 

9 8.5 and about 9.0; 

10 said classified category 4 corresponds with a predicted HbAu between about 

11 9.0 and about 9.5; 

12 said classified category S corresponds with a predicted HbAu between about 

13 9.5 and about 10.1; 

14 said classified category 6 corresponds with a predicted HbAu between about 

15 10.1 and about 11.0; and 

16 said classified category 7 corresponds with a predicted HbAu above about 

17 11.0. 

1 17. The system of claim 14, herein said processor being further progranuned to: 

2 define predicted confidence intervals for corresponding said HbAu , each of said 

3 predicted confidence intervals representing a range of values for HbAu* 
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1 18. The system of claim 1 7, v^eiein said predicted HbAic confidence intervals 

2 have about a 95% confidence level. 

1 1 9. A glycemic control system for evaluating HbAu of a patient, said system 

2 comprising: 

3 a BG acquisition mechanism, said acquisition mechanism configured to acquire BG 

4 data fiom the patient, 

5 a database component operative to maintain a database identifying said BG data; 

6 a processor progranomed to: 

7 compute weighted deviation toward high blood glucose (WR) and estimated 

8 rate of change of blood glucose (Dr) based on said collected BG data; and 

9 estimate HbAu using a predetermined mathematical formula based on said 

10 computed WR and Dr. 

1 20. A computer program product comprising a computer useable medium having 

2 computer program logic for enabling at least one processor in a computer system to evaluate 

3 HbAic of a patient based on BG data, said computer program logic comprising: 

4 computing weighted deviation toward high blood glucose (WR) and estimated rate of 

5 change of blood glucose (Dr) based on said collected BO data; and 

6 estimating HbAu using a predetermined mathematical formula based on said 

7 computed WR and Dr. 

1 21. The computer program product of claim 20, i^dierein said computer program 

2 logic further comprises: 

3 providing a predetermined confidence interval for classification of said estimated 

4 value of HbAiG, wherein said confidence interval is a single value or a range of values. 

1 22. A computerized method for evaluating the long term probability for severe 

2 hypoglycemia (SH) of a patient based on BG data collected over a predetermined duration, 

3 said method comprising: 
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4 computing weighted deviation toward low blood glucose (WL) and estimated rate of 

5 fall of blood glucose in the low BG range (DrDn) based on said collected BG data; and 

6 estimating the number of future SH episodes using a predetermined mathematical 

7 formula based on said computed WL and DrDn. 

1 23. The method of claim 22, wherein: 

2 said computed WL is mathematically defined fix>m a series of BG readings xu ... 

3 Xn taken at time points t2, ...» in as: 

4 »7:=iXw/(x,;2) 

5 where: 

6 wl(BG;a) = 10.f(BGf if f(BG)>0 and 0 otherwise, 

7 a =2, representing a weighting parameter, and 

8 said computed DR is mathematically defined as: 

9 DrDn « average of 5*+/ - 5*, provided that Sk< 

10 v^diere: ^lO.SOt-^ tif for *=0./. /« - //, 

11 S(tj)^f(Xj)j0Tj^l ft 

1 24. The method of claim 22, wherein said estimated number of future SH 

2 episodes (EstNSH) is mathematically defined as : 

3 EsfNSH = 3.3613(WL) - 4.3427(DrDn) - 1.2716. 

1 25. The method ofclaim 22, further comprising: 

2 defining predetenmned EstNSH categories, each of said EstNSH categories 

3 representing a range ofvalues for EstNSH; and 

4 assigning said EsfNSH to at least one of said EstNSH categories. 

5 26. The method of claim 25, herein said EstNSH categories are defined as 

6 follows: 

7 category 1, wherein said EstNSH category is less than about 0.775; 

8 category 2, wherein said EsfNSH category is between about 0.775 and about 3.750 ; 
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category 3, wherein said EstNSH category is between about 3.750 and about 7.000; 

and 

category 4, wherein said EstNSH category is above about 7.0. 

27. The method of claim 26, further comprising: 

defining a probability of incurring a select number of SH episodes respectively for 
each of said assigned EstNSH categories; wherein said probability and said respective select 
number of SH are defmed as: 

said classified category 1 corresponds with about a 90% probability of 
incurring about 0 SH episodes and about a 10% probability of mcurring about 1 or 
more SH episodes over the predetermined duration; 

said classified category 2 corresponds with about a S0% probability of 
incurring about 0 SH episodes, 25 % probability of incurring about 1 to about 2 SH 
episodes, and 25 % probability of incurring more than 2 SH episodes over the 
predetermined duration; 

said classified category 3 corresponds with about a 25% probability of 
incurring about 0 SH episodes, 25 % probability of incurring about 1 to about 2 SH 
episodes, and 50 % probability of incurring more than 2 SH episodes over the 
predetermined duration; and 

said classified category 4 corresponds with about a 20% probability of 
incurring about 0 to about 2 SH episodes and about a 80% probability of incurring 
more than 2 SH episodes over the predetermined duration. 

28. The method of claim 25, further comprising: 

defining a probability of incurring a select number of SH episodes respectively for 
each of said assigned EstNSH categories; and 

providing at least one probability of incuning a select number of SH episodes 
according to said EstNSH category to vAich said EstNSH is assigned. 
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1 29. A computerized method for evaluating the long term probability for severe 

2 hypoglycemia (SH) of a patient based on BG data collected over a predetermined duration, 

3 said method comprising: 

4 computing weighted deviation toward low blood glucose (WL) and estimated rate of 

5 fall of blood glucose in the low BG range (DrDn) based on said collected BG data; 

6 estimating the number of future SH episodes using a predetermined mathematical 

7 formula based on said computed WL and DrDn; and 

8 defining a probability of incurring a select number of SH'episodes respective to said 

9 estimated SH episodes. 

1 30. A system for evaluating the long term probability for severe hypoglycemia 

2 (SH) of a patient based on BG data collected over a predetermined duration, said system 

3 comprising: 

4 a database component operative to maintain a database identifying said BG data; 

5 a processor programmed to: 

6 computing weighted deviation toward low blood glucose (WL) and estimated 

7 rate of fall of blood glucose in the low BG range (DrDn) based on said collected BG 

8 data; and 

9 estimating the number of future SH episodes using a predetermined 

10 mathematical formula based on said computed WL and DrDn. 

1 31. The system of claim 30, wherein: 

2 said computed WL is mathematically defined from a series of BG readings ... 

3 x„ taken at time points 0, tj, /« as: 

4 WL-^f^wKxra) 

5 where: 

6 wl(BG;a) = 10f(BGf if f(BG)>0 and 0 otherwise, 

7 a »2, representing a weighting parameter, and 
8 
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9 said computed DR is mathemadcally defined as: 

10 DrDn = average of j*+7 - su, provided that Sk< 5*+/, 

1 1 where: 

12 Sk=10.S(k-^ti)^ fox k^0.l,..,U-th 

13 S(tj)^f(xj\totj^l...,n. 

1 32. The system of claim 30, viierein said estimated number of future SH episodes 

2 (EstNSH) is mathematically defined as : 

3 EstNSH = 3.3613(WL) - 4.3427(DrDn) - 1.2716. 

1 33. The system of claim 30»A^erein said processor being furdier programmed to: 

2 define predetermmed Esfl^SH categories, each of said EstNSH categories 

3 representing a range of values for EstNSH; and 

4 assign said EstNSH to at least one of said EstNSH categories. 

1 34. The system of claim 33, >;^rein said EstNSH categories are defined as 

2 follows: 

3 category 1 , wherein said EstNSH category is less than about 0.77S; 

4 category 2, wherein said EstNSH category is between about 0.77S and about 3.750 ; 

5 category 3, wherein said EstNSH category is between about 3.7S0 and about 7.000; 

6 and 

7 category 4, wdierem said EstNSH category is above about 7.0. 

1 35. The method of claim 34, wherein said processor being fiirther programmed to: 

2 define a probability of incurring a select number of SH episodes respectively for each 

3 of said assigned EstNSH categories, wherein said probability and said respective select 

4 number of SH are defined as: 

5 said classified category 1 corresponds with about a 90% probability of 

6 incurring about 0 SH episodes and about a 10% probability of incurring about 1 or 

7 more SH episodes over the predetermined duration; 
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8 said classified category 2 corresponds with about a 50% probability of 

9 incurring about 0 SH episodes, 25 % probability of incurring about 1 to about 2 SH 

10 episodes, and 25 % probability of incurring more than 2 SH episodes over the 

1 1 predetermined duration; 

12 said classified category 3 corresponds with about a 25% probability of 

13 incurring about 0 SH episodes, 25 % probability of incuning about 1 to about 2 SH 

14 episodes, and 50 % probability of incurring more than 2 SH episodes over the 

1 5 predetermined duration; and 

16 said classified category 4 corresponds with about a 20% probability of 

17 incurring about 0 to about 2 SH episodes and about a 80% probability of incurring 

1 8 more than 2 SH episodes over the predetermined duration. 

1 36. The system of claim 33, wherein said processor being further programmed to: 

2 define a probability of incurring a select number of SH episodes respectively for each 

3 of said assigned EstNSH categories; and 

4 provide at least one probability of incurring a select number of SH episodes according 

5 to said EstNSH category to which said EsfNSH is assigned. 

1 37. A glycemic control system for evaluating the long term probability for severe 

2 hypoglycemia (SH) of a patient, said system comprising: 

3 a BG acquisition mechanism, said acquisition mechanism configured to acquire BG 

4 data from the patient, 

5 a database component operative to maintain a database identifying said BG data; 

6 a processor programmed to: 

7 computing weighted deviation toward low blood glucose (WL) and estimated 

8 rate of ML of blood glucose in the low BG range (DrDn) based on said collected BG 

9 data; and 

10 estimating the number of future SH episodes using a predetermined 

1 1 mathematical formula based on said computed WL and DrDn. 
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1 38. A computer program product comprising a computer useable medium having 

2 computer program logic for enabling at least one processor in a computer system to evaluate 

3 long term probability for severe hypoglycemia (SH) of a patient based on BG data, said 

4 computer program logic comprising: 

5 computing weighted deviation toward low blood glucose (WL) and estimated rate of 

6 fall of blood glucose in the low BG range (DrDn) based on said collected BG data; and 

7 estimating the number of future SH episodes using a predetermined mathratiatical 

8 formula based on said computed WL and DrDn. 

1 39. The computer program product of claim 38, herein said computer program 

2 logic further comprises: ( 

3 defining a probability of incurring a select number of SH episodes respective to said 

4 estimated SH episodes. 

1 40. A computerized method for evaluating the short term risk for severe 

2 hypoglycemia (SH) of a patient based on BG data collected over a predetermined duration, 

3 said method comprising: 

4 computing weighted deviation toward low blood glucose (WL); 

5 detennining Maxfwl) by calculating maximum value of wl(BG;2); and 

6 determining risk value by taking the geometric mean of WL and Max(wl) over said 

7 predetermined duration, said risk value is mathematically defined as: 

8 riskvalue - ^f¥L9Max{wl) . 

1 41 . The method of claim 40, wherein: 

2 said computed WL is mathematically defined from a series of BG readings xi, xt, ... Xn 

3 taken over the predetermined duration as: 

4 WL-^f^wl{x,:i) 

5 where: 

6 wl(BG:a) = lO/fBG)" if f(BG)>0 and 0 otherwise, 

7 a -2, representing a weighting parameter. 
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1 42. The method of claim 40, further comprising: 

2 providing a predetermined threshold risk value; and 

3 comparing said determined risk value to said threshold risk value. 

1 43 . The method of claim 42, wherein: 

2 if said determined risk value is greater than said threshold value then short term risk 

3 of incurring a hypoglycemic episode is high; and 

4 if said determined risk value is less than said threshold value then short term risk of 
s mcurring a hypoglycemic episode is low. 

1 44. The method of claim 43, herein said short term is approximately a 24 hoxir 

2 period. 

1 4S. The method of claim 43, wherein said short term ranges from about 12 to 

2 about 72 hour period. 

I 46. The method of claim 43, herein said threshold value is approximately 17. 

1 47. The method of claun 43, wherein said threshold value is between about 12 to 

2 25. 

1 48. A system for evaluating the short term risk for severe hypoglycemia (SH) of a 

2 patient based on BG data collected over a predetermined duration, said system comprising: 

3 a database component operative to maintain a database identifying said BG data; 

4 a processor programmed to: 

5 compute weighted deviation toward low blood glucose (WL); 

6 determine Max(wl) by calculating maximum value of wl(BG;2)\ and 

7 determine risk value by taking the geometric mean of WL and Max(wl) over 

8 said predetermined duratioa said risk value is mathematically defined as: 

9 risk value = ^WL • Max{wl) . 
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1 49. The system of claim 48, wherein: 

2 said computed WL is mathemAticaily defined from a series of BG readings xu ... 

3 taken over the predetermined duration as: 

5 where: 

6 wl(BG;a) = 10.f(BG/ if f(BG)>0 and 0 otherwise, 

7 a »2, representing a weighting parameter. 

1 SO. The system of claim 48, wherein said processor being further programmed to: 

2 provide a predetermined threshold risk value; and 

3 compare said determined risk value to said threshold risk value. 

1 SI. The system of claim SO, wherein: 

2 if said determined risk value is greater than said threshold value then short term risk 

3 of incurring a hypoglycemic episode is high; and 

4 if said determined risk value is less than said threshold value then short term risk of 
s incurring a hypoglycemic episode is low. 

1 S2. The system of claim SI, viierein said short term is approximately a 24 hour 

2 period. 

1 S3. The method of claim SI, wherein, said short term ranges from about 12 to 

2 about 72 hour period. 

1 S4. The system of claim S 1 , wherein said threshold value is approximately 1 7. 

1 SS. The method of claim SI, wherein said threshold value is between about 12 to 

2 2S. 
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1 56. A glycemic control system for evaluating the short term risk for severe 

2 hypoglycemia (SH) of a patient, said system comprising: 

3 a BG acquisition mechanism, said acquisition mechanism configured to acquire BG 

4 data from the patient, 

5 a database component operative to maintain a database identifying said BG data; 

6 a processor programmed to: 

7 compute weighted deviation toward low blood glucose (WL); 

8 determine Max(wl) by calculating maximum value of wl(BG;2); and 

9 determine risk value by taking the geometric mean of WL and Max(wl) over 

10 said predetermined duratioa said risk value is mathematically defined as: 

11 risk value = ^WL • Mbc(w/) . 

1 57. A computer program product comprising a compiiter useable medium having 

2 computer program logic for enabling at least one processor in a computer system to evaluate 

3 the short term risk for severe hypoglycemia (SH) of a patient based on BG data collected over 

4 a predetermined duration, said computer program logic comprising: 

5 computing weighted deviation toward low blood glucose (WL); 

6 determining Maxfwl) by calculating maximum value of wl(BG;2); and 

7 determining risk value by takmg the geometric mean of WL and Max(wl) over said 

8 predetermined duration, said risk value is mathematically defined as: 

9 risk value - ^IWL^Max^wl) , 

1 58. The computer program product ofclaim 57, wherein said computer program 

2 logic fiirther comprises: 

3 providing a predetermined threshold risk value; and 

4 comparing said determined risk value to said threshold risk value, 
s 
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